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Abstract. We present results of a survey of the spectra 
of Low Mass X-ray Binaries using ASCA. It is shown that 
all sources in the survey are well-fitted by the same two- 
component emission model that we have previously shown 
is able to describe both the non-dip and dip spectra of 
the dipping class of LMXB. This model consists of point- 
like blackbody emission from the neutron star plus Comp- 
tonized emission from a disk-like accretion disk corona of 
radius typically 50,000 km. Additional data from results 
published elsewhere by us from BeppoSAX and ASCA are 
added to the survey. The large variation in blackbody lu- 
minosity of survey sources is shown to be due primar- 
ily to major changes in blackbody emitting area. Fitting 
a multi-temperature disk blackbody plus Comptonization 
model to the survey spectra requires values of inner disk 
radius substantially less than the neutron star radius in 
many cases, making disk origin of the blackbody highly 
unlikely. Assuming that the emission is from an equatorial 
strip on the neutron star, it is shown that the half-height 
of the strip h agrees well with the half-height H of the 
radiatively-supported inner accretion disk, this agreement 
spanning three orders of magnitude in each parameter. 
Possible mechanisms for the agreement are discussed, in- 
cluding radial accretion flow between inner disk and star, 
and accretion flow "creep" on the surface of the neutron 
star. 
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1. Introduction 

Low mass X-ray binaries (LMXB) may be classified ac- 
cording to inclination angle and according to colour-colour 
properties. Frank et al. (1987) explained the markedly 
different light curves of LMXB in terms of a model of 
varying inclination. Sources having inclination angles be- 
tween zero and ~ 65 will show no orbital related vari- 
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ability, sources with inclination of 65-80 will exhibit X- 
ray dipping, and sources with inclinations of 80-90 will 
be Accretion Disk Corona (ADC) sources. Strong orbital- 
related variability is seen in the light curves of the dip- 
ping and ADC sources, but not in the other sources which 
are viewed from above the orbital plane. The second type 
of classification is based on examining the behaviour of 
LMXB on X-ray colour-colour diagrams. Hasinger et al. 
(1989) first used the term "Z-track" to describe the evo- 
lution on a colour-colour diagram of bright LMXB such 
as ScoX-1 and CygX-2, the sources clearly undergoing 
marked spectral changes as they moved along the hor- 
izontal, normal and flaring branches. The Atoll sources 
were similarly recognised as having different colour-colour 
tracks (Hasinger & van der Klis 1989). Associated with 
tracks in the Z-track and Atoll sources was correlated tim- 
ing behaviour, especially of QPOs (e.g. van der Klis 1995). 

Clearly, strong spectral changes take place during track 
movement, but these are not understood. A number of sur- 
veys of the spectra of the Atoll and Z-track sources have 
been made, including that based on Exosat data (White 
et al. 1988). Schulz et al. (1989) classified a large number 
of sources in terms of colour-colour diagrams, and carried 
out spectral fitting using a blackbody plus cut-off power 
law model. A similar blackbody plus power law model was 
used by Schulz (1999) in a Rosat survey of LMXB. Chris- 
tian & Swank (1997) produced a compilation of the prop- 
erties of LMXB from the Einstein solid state and monitor 
proportional counter detectors. 

There are considerable differences between the emis- 
sion models that have been proposed for LMXB in general. 
A disk blackbody plus a Comptonized blackbody model 
was proposed by Mitsuda et al. (1989). A model also 
based on accretion disk blackbody emission was applied 
by Czerny et al. (1986) to Einstein data. The generalised 
thermal model (White et al. 1985) was used by White 
et al. (1988), with a blackbody added for high-luminosity 
sources which was not required for the lower luminosity 
burst sources, from which they concluded that spectral 
formation was dominated by Comptonizaton in the inner 
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accretion disk. More recently, a two-component model pro- 
posed by Church and Baluciiiska-Church (1995) has been 
shown to explain spectral evolution in the dipping LMXB 
sources, consisting of point-source blackbody emission 
identified with the neutron star plus Comptonized emis- 
sion from an extended ADC modelled by a cut-off power 
law (Church & Bahicinska-Church 1993, 1995; Church et 
al. 1997, 1998a,b; Balucinska- Church et al. 1999, 2000, 
2001). 

For the dipping LMXB sources, spectral analysis is 
more strongly constrained because of the requirement that 
a model fits not only the non-dip spectrum but also sev- 
eral levels of dipping. A major advantage of investigat- 
ing spectral evolution in these sources is that dipping is 
characterised by the slow removal of Comptonized emis- 
sion, and measurement of dip ingress times proves that the 
Comptonizing ADC has radial extent typically 50,000 km 
(Church 2000), which is strongly inconsistent with Comp- 
tonization models involving a small region in the locality 
of the neutron star, allowing such models to be discounted. 

Table 1. The observations 



Source 


Type 


Date 




GX9-h9 


Atoll 


1994, 


Mar 25 


GX13+1 


Atoll 


1994, 


Sep 10 


4U 1636-536 


Atoll 


1993, 


Aug 9 


SerX-1 


Prob Atoll 


1994, 


Sep 20 


X2127-hll9 


Prob Atoll 


1995, 


May 16 


AqIX-1 




1994, 


Apr 30 


X 1746-371 


Dipping 


1995, 


Sep 21 


XB 1254-690 


Dipping 


1994, 


Mar 18 


GX5-1 


Z-track 


1993, 


Oct 9 


CygX-2 


Z-track 


1993, 


Jun 10 



The success of this emission model with the dip- 
ping sources suggests that the model be tested with the 
other classes of LMXB, i.e. the Z-track and Atoll sources 
specifically. Other authors have successfully applied two- 
component blackbody -|- Comptonization models to Z- 
track and Atoll sources (above, and e.g. Hasinger et al. 
1990, Barret et al. 2000). However, these often ascribe 
emission regions to the blackbody and Comptonized emis- 
sion such as the accretion disk and a localised inner region 
respectively, e.g. Barret et al. (2000), and so physically, 
these models differ radically from ours and are also incon- 
sistent with our measured ADC sizes. Recently, the com,ptt 
implementation of the Titarchuk (1994) Comptonization 
model has been used (together with a blackbody term) 
to fit BeppoSAX data on LMXB sources (e.g. in't Zand 
et al. 1999). It has been suggested that comptt is prefer- 
able to the cut-off power law which may overestimate the 
spectrum at low energies because of lack of soft seed pho- 
tons. However, it is easy to calculate the spectrum of the 



disk (expressed as the photon number flux per kcV) inte- 
grated to a radius equal to our measured ADC sizes, by 
combining the temperature gradient T(r) from standard 
thin disk theory with Planck's equation. This shows that 
for typical LMXB luminosities, the integrated spectrum 
forms a broad peak between 0.001 and 0.1 keV, thus pro- 
ducing a huge sea of soft photons, so that the cut-off power 
law is perfectly good above ~0.1 keV. However, fitting the 
comptt model has produced high values of average kT for 
the seed photons, e.g. 1.0-2.0 keV (Guainazzi et al. 1998) 
which are >> than the above expected values and incon- 
sistent with the assumption of Titarchuk (1994) that the 
seed photons are very soft: kT « kTc (the electron tem- 
perature). This will lead to a substantial underestimation 
of the Comptonized spectrum at ^1 keV as the Wien ap- 
proximation (good for hv » kT) is used in comptt to 
describe the seed photon spectrum. 

In this paper, we apply the blackbody plus cut-off 
power law model to as many non-dipping Atoll and Z- 
track sources observed using ASGA as possible plus the 
non-dip spectra of two dipping sources. We show that the 
model provides a good description of the spectra of these 
sources. 

2. Analysis and Results 

The observations of 10 LMXB with ASCA (Tanaka et al. 

1994) were made between 1993, June and 1995, September 
as shown in Table 1. The designation "probably atoll" 
was taken from van der Klis (1995) when the type is not 
well known. For each source, raw data files containing less 
than 100 counts were deleted, and the temperature-based 
gain correction was made. Screening was carried out after 
inspection of plots of the housekeeping parameters. Data 
were selected with the source less than 0.01 from the 
telescope pointing direction, elevation above the Earth's 
rim more than 5 and geomagnetic rigidity more than 6 
GeV/c. Particle background was rejected for values of the 
background monitor rate of more than 200 count s~^ in 
most cases. Source data were selected from a 6' radius 
circular region in the image. Many of the sources were 
bright and so background subtraction was not important; 
when necessary, background data were extracted from an 
appropriate 6' region. Deadtime correction was performed 
on both light curves and spectra; this enabled high and 
medium bit rate data to be combined in a single light 
curve. In all cases except one, GIS2 data were used; in 
the case of X 1746-371, GIS3 data were used because of 
problems with the GIS2 raw data. 

Light curves were extracted in the total CIS energy 
band and also in 3 standard bands: 0.7-4.0 keV, 4.0-7.0 
keV and 7.0-10.0 keV. These were used to examine source 
variability and to construct colour-colour diagrams using 
the ratio of count rates 7.0-10.0 keV/4.0-7.0 keV and 4.0- 
7.0 keV/0. 7-4.0 keV. Generally, there was little evidence 
for movement on the colour-colour diagrams during obser- 
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Fig. 1. Best fits to the ASCA spectra of 10 LMXB using the two-component continuum model discussed in the text, 
either in the form of blackbody + power law or blackbody + cut-off power law. Lines were added to the model in the 
4 sources where definite line detections were made. In all cases the ordinate is in units of photon cm~^ s~^ keV~^ 



vations; some of the soiirccs were constant, some showed 
slow drifts in intensity and some more rapid variability. 
To ensure that data could not be mixed from different 
intensity states, spectra were selected from appropriately 
narrow intensity bands for all sources, for example, in the 
band 75-80 count s~^ for GX13+1. Response files were 
constructed appropriate to the position of each source in 



the image, and the latest corrections applied to the re- 
sponses to correct for known differences between the de- 
tectors. 

Spectra were grouped to a minimum of 100 counts per 
bin for strong sources and to 20 counts per bin for faint 

sources, and systematic errors of 2% added (Fukazawa 
et al. 1997). For each source, a range of models was 
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Table 2. Spectral fitting results. In most cases, the blackbody + cut-off power law model did not give any improvement 
over the blackbody + power law model, indicating a cut-off energy much higher than 10 keV. Only in cases where 
there was an improvement are results for both of these models given. 



Source 


Model 


iVu 
10^2 cm-2 


kT 
kev 


r 


keV 


xVdof 


Atoll sources 


GX9+9 


Power law 


0.55 ±0.01 




1.89 ±0.07 




4171/712 




Cut-off power law 
Blackbody -|- power law 
Blackbody + cut-off power law 


0.19 ±0.01 
0.37 ±0.02 
0.26 ±0.01 


1.17 ±0.02 
1.27 ±0.06 


0.80 ±0.03 
1.95 ±0.03 
1.21 ±0.03 


4.1 ±0.1 
~6 


507/711 
615/710 
514/710 


GX13-f-l 


Power law 
Cut-off power law 
Blackbody -I- power law 


4.4 ±0.1 

2.4 ±0.1 

3.5 ±0.2 


1.18 ±0.02 


2.51 ±0.02 
-0.30 ±0.13 

2.88 ±0.22 


1.9 ±0.1 


1639/502 
459/501 

446/499 


4U 1636-536 


Power law 
Cut-off power law 
Blackbody -|- power law 


0.55 ±0.02 
0.46 ± 0.03 
0.43 ± 0.05 




2.12 ±0.02 
1.83 ±0.09 
1.99 ± 0.07 


13.511^ 


423/396 
395/395 
391/394 


SerX-1 


Power law 

Cut-off power law 
Blackbody -|- power law 


0.71 ± 0.02 

0.43 ±0.02 
0.62 ± 0.03 


1.26 ± 0.08 


1.85 ±0.02 
1.05 ±0.05 
2.00 ±0.05 


5.4 ±0.4 


1253/605 
513/604 
514/603 


X 2127+119 


Power law 

Cut-off power law 
Blackbody -I- power law 


0.58 ± 0.04 

0.14 ±0.06 
0.24 ±0.09 


0.89 ± 0.06 


2.04 ± 0.04 

0.56 ±0.19 
1.85 ±0.16 


9 7+0.4 
^•'-0.3 


400/211 

218/210 
216/209 


Aql X-1 


Power law 
Cut-off power law 
Blackbody -|- power law 


0.64 ± 0.03 
0.48 ± 0.05 
0.49 ± 0.08 


n QO-l-O-25 
"-"■^^-0.10 


1.75 ±0.02 
1.39 ±0.10 
1.66 ± 0.08 




535/481 
497/480 
503/479 


Dipping sources 


XB 1746-371 


Power law 


0.50 ± 0.05 




1.83 ±0.06 




129/106 




Cut-off power law 
Blackbody -|- power law 


0.31 ±0.09 
0.50 ±0.05 


-1 99+0.31 
i-^^_0.21 


1.20 ±0.27 
1.91 ±0.05 


"•^-2.0 


115/105 
283/276 


XB 1254-690 


Power law 
Cut-off power law 
Blackbody -|- power law 


0.48 ±0.03 
0.21 ± 0.05 
0.38 ± 0.06 


1.17±0.11 


1.83 ±0.03 
0.93 ±0.13 
1.96 ±0.13 


^•"-0.6 


389/281 
259/278 
264/277 


Z-track sources 


GX5-1 


Power-law 


4.8 ±0.1 




2.56 ± 0.03 




952/438 




Cut-off power law 
Blackbody -|- power law 


3.0 ±0.1 
3.4 ±0.2 


1.05 ±0.04 


0.20 ±0.18 
2.29 ±0.18 


2.2 ±0.2 


476/437 
475/436 


CygX-2 


Power- law 
Cut-off power law 
Blackbody -|- power law 
Blackbody -|- cut-off power law 


0.56 ±0.01 
0.06 ± 0.02 
0.29 ± 0.02 
0.10 ±0.02 


1.01 ±0.02 
1.06 ±0.12 


2.11 ±0.01 
0.45 ± 0.05 
2.16 ±0.05 
0.75 ± 0.06 


2.5 ±0.1 
~3 


3186/541 
483/540 
559/539 
487/539 



tried including simple one-component models such as 
absorbed power law (ab*PL, where AB is the absorp- 
tion term and PL the power law), absorbed blackbody 
(ab*bb) and absorbed bremsstrahlung (ab*br). The ab- 
sorbed blackbody model always gave totally unacceptable 
fits with x^/degree of freedom (dof) typically 10; for the 
bremsstrahlung models the fits were also not acceptable 
with %^ /dof varying between 1.2 and 2.0. Two-component 
models consisting of blackbody plus a Comptonized com- 
ponent were tried, in the forms ab*(bb ± pl) and ab*(bb 
+ cpl). The first form was used to approximate a cut- 



off power law at energies well below the Comptonization 
break. In several cases discussed below, it was necessary to 
add lines to the continuum model to obtain an acceptable 
fit. Results for the absorbed power law, absorbed cut-off 
power law and two-c;omponent models are given in Table 
2 and these were obtained with lines added if appropriate. 

In most sources, there was clear curvature, i.e. down- 
curving, in the spectrum between 1 and 10 keV (see 
Fig. 1), such that the ab*pl model gave poor fits as 
can be seen from the values in Table 2. The gen- 
eralised thermal (model AB*CPL, where CPL is a cut-off 
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Table 3. Line detections 



Source Energy Width a EW 

keV keV eV 



GX13+1 


6.37 ±0.23 


0.1 


44 


SerX-1 


6.6 


0.17 


81 


AqlX-1 


0.99 ± 0.04 


0.01 


53 




6.81 ±0.16 


0.001 


43 


XB 1254-690 


6.40 ±0.12 


0.001 


98 


law), provided 


good fits to 


all of the 


spectra. 



ever, in the cases of several sources, particularly GX 13+1, 
X 2127+119, GX5-1 and CygX-2, this was only achieved 
with very small and even negative values of the power law 
index T which may be regarded as improbable. The values 
of cut-off energy Eqo are all very small (^few keV) which 
is typical of this model. However, we know from work on 
several dipping sources using BeppoSAX that the break 
energy is often much greater than 10 keV (Church et al. 
1998b; Bahicinska-Church et al. 1999, 2000) and so the 
low values in most cases are unlikely to be real, but result 
from this model fitting the curvature between 1 and 10 
keV actually due to the presence of a blackbody. Fitting a 
two-component model showed that in all cases, either the 
blackbody plus power law model or the blackbody plus 
cut-off power law model gave a fit as good or better than 
the AB*CPL model, but without requiring any odd values. 
In most cases, the simpler ab*(bb + PL) model gave a 
fit equally as good as the AB*(bb + CPl) model showing 
that Eco was >> 10 keV and so not possible to be deter- 
mined using ASCA. In the case of GX9+9 and CygX-2 
this was not the case, showing that Eco was less than 10 
keV. Fitting for these sources, was carried out with Eqo 
fixed at a sequence of values as discussed below. 

Lines were detected in the sources GX13+1, AqlX-1, 
SerX-1 and XB 1254-690, and the energies, widths and 
equivalent widths are shown in Table 3. These agree gen- 
erally with the detections made in the extensive study of 
lines in LMXB using ASCA by Asai et al. (2000). There 
was also evidence for a weak line in XB 1746-371 at ~6.4 
keV also detected in SIS by Asai et al. , and a possible weak 
detection of a line at ~1.5 keV in GX 13+1. In SerX-1, the 
iron feature detected was broad, and the width was not 
well-constrained using GIS data; consequently, to obtain 
the best continuum fit, the line energy and width a were 
set to the values determined by Asai et al. The best-fit 
two-component model for each source is shown in Fig. 1. 

2.1. X 1746-371, XB 1254-690 

Spectral analysis results of the ASCA observations of these 
dipping sources have not previously been presented. These 

sources were well fitted by the two-component model. In 
XB 1746-371, the best- fit results were obtained by simul- 



taneous fitting of the non-dip spectrum and the deepest 
dip spectrum. This technique has previously been used to 
constrain emission parameters better in dipping sources 
(e.g. Balucinska-Church et al. 1999). 

2.2. GX9+9, Cyg X-2 

In GX9+9 and Cyg X-2 the blackbody plus power law 
model gave a substantially worse fit than the single com- 
ponent cut-off power law model showing that Eqo is not 
>> 10 keV. The model AB*(bb + CPl) gave improved 
qualities of fit. It was difficult to determine the cut-off en- 
ergy because of the competition in modelling the 1-10 keV 
curvature between the blackbody and the Comptonization 
down-curving, and so Eco was fixed at a series of values: 
10, 8, 6, 4 and 2 keV. In GX9+9, became as good 
as for the cut-off power law model for Eqo ~6 keV. For 
smaller cut-off energies, continued to improve, however 
the power law index decreased to < 1 which is unlikely to 
be real. For Cyg X-2, stepping in Eqo gave acceptable 
for Eco ~3 keV. This cut-off energy agrees well with the 
value of 3.32 keV obtained by Smale et al. (1993) from 
BBXRT. Thus, our preferred values of cut-off energy are 
~6 keV and ^3 keV in GX9+9 and Cyg X-2 respectively. 
Clearly, these should be re-measured in a wider band. 
However, test showed that the results presented below are 
not sensitive to the value of Eco- 

2.3. Disk blackbody modelling 

Finally, we tested the model consisting of disk blackbody 
emission plus a Comptonization term in the form of a 
cut-off power law, adding lines as before to the sources re- 
quiring lines. The 10 sources studied in the present ASCA 
survey were fitted, plus 3 sources recently studied using 
BeppoSAX (Balucinska-Church et al. 1999, 2000; Church 
et al. 1998b). Results shown in Table 4 include values of 
inner disk radius obtained from the normalization of 
the blackbody via: r\ = (norm /cos i)^/^ • rf, where i is the 
inclination angle of the binary system and d is the source 
distance. Inclination angles were assumed of - 70 for 
Atoll and Z-track sources, 70 - 80 for dipping sources, 
except for X 1624-490 where 60 65 was assumed, the 
upper limit being set by the absence of X-ray eclipses. 
For Cyg X-2, we use the range 58.5 - 66.5 (Orosz & Ku- 
ulkers 1999). Thus for each source a range of n values 
is shown, with 90% confidence uncertainties derived from 
the normalization uncertainties shown for the larger value 
in each case. Although good or acceptable fits were ob- 
tained, the results in many cases had improbably small or 
negative values of power law index. The values of ri were 
in 8 out of 13 cases unphysical, i.e. < the radius of the neu- 
tron star i?*, in 5 of these, n being 10 times smaller than 
R^. In 5 other cases, GX9+9, GX13+1, GX5-1, Cyg X-2 
and 4U 1636-536, n was not clearly unphysical, but the 
power law index was small or negative. We also note a 
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Table 4. Results of fitting a disk blackbody + cut-off power law model to ASCA and BeppoSAX data. The predomi- 
nance of small or negative values of photon index T can be seen, and the unphysical values of n less than the radius 
of the neutron star. 



Source 


Instrument 


keV 


r 


Eco 
keV 


xVdof 


r\ 

km 


GX9-h9 


ASCA CIS 


0.70±0.10 


-ooi+"-^'' 


9 7+(l.a 
^- ' -0.4 


490/709 


10-^17+^ 


GX 13-M 


ASCA CIS 


1 14+0.32 
■'^•J^^-0.23 


-3.00±0.01 


1.17±0.01 


435/499 


ll-19tj° 


4U 1636-536 


ASCA GIS 


0.64±0.07 


+0.50lgi^ 


4.0±0.8 


365/393 


14-24t^ 


ScrX-1 


ASCA GIS 


2 42+"-i=^ 

^•^^-0.08 


+0.20«;«J 




489/601 


2.7-4.7t°-j 


X 2127+119 


ASCA GIS 


3.22±1.24 




1 5+°-2 


213/208 


0.79-1. l+![? 


AqlX-1 


ASCA GIS 


1.00±0.21 


-1.05±1.62 


9 Q+3.3 


482/480 


2.8-4.8tJ-j' 


XB 1746-371 


ASCA GIS 


q no+S.Oe 


.1 14+0.58 
-'-•-'-^-1.67 


0.85±0.01 


112/103 


1.0-1.4lJ;f 


XB 1254-690 


ASCA GIS 


9 CI +0.60 
"-^ 1 93 


-hl.42±0.63 


3.98±0.1 


257/276 


2.0-2.9±2.1 


GX 5-1 


ASCA GIS 




-1.36±0.10 




473/435 


15-26t2 


CygX-2 


ASCA GIS 


1.83±0.17 


-0.90±0.04 


0.74±0.01 


480/538 


8.6-9.813.0 


XB 1916-053 


SAX NFI 


2.96±0.24 


+1.60±0.04 


200±150 


556/520 


0.48 0.67±0.13 


XB 1323-619 


SAX NFI 


2 04+0. 56 


.7+0.24 

T^J--J- ' -0.39 


32±11 


204/216 


0.48-0.67±0.01 


X 1624-490 


SAX NFI 


5 1+1-6 


-0.78±0.58 


1.7±0.3 


291/294 


0.39-0.4lt°;|5 



tendency using ^S'CVl to underestimate kTj and overesti- 
mate ri compared with BeppoSAX as shown by restricting 
the BeppoSAX band to 1-10 keV, so in some cases, the ri 
values in Table 4 may be too large. It is known that the 
multi-temperature disk blackbody model used can under- 
estimate Ti, however, maximum correction factors of ~2-4 
may be calculated based on the published work of Merloni 
et al. (2000) and Kubota et al. (2000). Allowing for this, 
the very small values of ri in Table 4 make it unlikely that 
the blackbody emission is from the disk. 

We have shown that the two-component model (either 
as BB -|- PL or as bb -|- CPl) gives very good fits to all of 
the sources in our survey. In all blackbody com- 

ponent is required, and we next examine the properties of 
this component. 

2.4^. The Blackbody- Total Luminosity Relation 

In Fig. 2 the variation of blackbody luminosity Lbb 
in the band 1-30 keV with the total luminosity 
Z/Tot in the same band is shown. Also plotted in Fig. 
2 are additional points from our analyses of the 5ep- 
poS'.dX observations of XB 1916-053 (Church et al. 1998b), 
XB 1323-619 (Balucihska-Church et al. 1999) and X 1624- 
480 (Balucihska-Church et al. 2000) and for XBT 0748-676 
from ASCA data (Church et al. 1998a). These points are 
useful because the sources arc weaker, and spectral pa- 
rameters well-determined from the broadband BeppoSAX 
data. Source distances were taken from the Einstein Ob- 
servatory survey of LMXB of Christian & Swank (1997). 
In the case of XB 1323-619, a distance of 10 - 20 kpc 
was derived by Parmar et al. (1989) based on the mea- 
sured peak burst flux assumed to be Eddington-limited. 
As this may have underestimated the peak flux, we have 



assumed a value of 10 kpc. We also assumed 10 kpc for 
XBT 0748-676, the distance of which is not known. The 
sources in the present survey cover the luminosity range 
^ 5 X 10"^® to 5 X 10"^^ erg s^^. Errors in Xbb arc derived 
from the 90% confidence limits in the blackbody normal- 
ization; similarly the errors in iTot were obtained from the 
90% errors in power law normalization. Possible errors in 
source distance will move points equal amounts on both 
axes. 

Additionally, a line is drawn showing the Newtonian 
value of 50% of Lxot that should remain at the inner disk 
and be available for neutron star emission. In fact, the 
ratio of the energy available to the star to the energy in 
the disk depends on the spin of the neutron star (Sunyaev 
& Shakura 1986), varying from 2.2 for a slowly spinning 
star to small values for a star spinning rapidly, implying in 
the first case that the star can receive 69% of the energy. 
It can be seen that there is systematic behaviour in Fig. 
2 such that the brighter sources approach the 50% line 
whereas the weaker sources fall increasingly below this 
line. Schulz et al. (1989) and Hasinger et al. (1990) found 
a similar variation of the blackbody luminosity with the 
total luminosity in CygX-2. It should be stressed that 
we do not claim a simple relationship between Xbb and 
Lxot, but rather a broad band in the Figure where the 
sources lie. We next examine the reasons for the variation 
in blackbody luminosity over the sources investigated. 

2.5. Blackbody Properties 

Table 5 lists the blackbody parameters derived from spec- 
tral fitting. Included here are the blackbody temperature 
fcTeB, the blackbody luminosity in the band 1-30 keV 
and the blackbody radius i?BB defined by the relation Xbb 
= A-K R^^aT^, where a is Stefan's constant. Also shown 
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XB 1354-690 



XBT 0748-676 



T4U 1636-536 



XB 1916-053 TXB 1746-371 



XB 1323-619 



~|~Aql X-1 



10 100 

Lum^^j (1 30 kcV) 



Fig. 2. Variation of blackbody luminosity in the band 1-30 keV with total luminosity in the same band; both lu- 
minosities are plotted in units of 10^^ erg s~^. The dotted line shows the relationship expected from simple energy 
considerations that the blackbody luminosity of the neutron star should equal the luminosity of the accretion disk; 
i.e. 50% of the total luminosity 



is the half-height h of the emission region assumed to be 
an equatorial strip on the neutron star. It can be seen 
that Lbb is reasonably well correlated with Rbb- In Fig. 
3 we show the variation of blackbody temperature with 
blackbody radius. It is striking that the major cause of 
variation is the emitting area. As the Lxot increases as we 
move from the the extreme left of Fig. 3 (XB 1323-619) 
to the extreme right (GX5-1), kT^s falls by 40%, so that 
decreases by a factor of 8. Rbb however, increases by 
a factor of ~75 and emitting area by a factor of 5400 so 
that Lbb increases by a factor of ^700, consistent with 
the change of Lbb between these sources in Fig. 2. 



Table 5. Blackbody properties, h is the half- height of 
the emission region on the neutron star; the blackbody 
luminosity Lgg (1-30 keV) is given in units of 10^^ erg 



Source 


d 


kTBB 


7-36 

^BB 


-Rbb 


h 




kpc 


keV 




km 


km 


GX9+9 


5.0 


1.27 


1.66 


3.2 


1.03 


GX 13+1 


7.0 


1.18 


39.5 


12.7 


16.0 


4U 1636-536 


6.5 


0.65 


1.52 


8.6 


7.5 


SerX-1 


8.4 


1.26 


15.7 


7.0 


4.9 


X 2127-^119 


13.0 


0.89 


3.67 


7.0 


4.9 


AqlX-1 


4.8 


0.92 


0.29 


1.81 


0.33 


XB 1746-371 


9.0 


1.22 


0.58 


1.4 


0.21 


XB 1254-690 


12.0 


1.17 


4.74 


4.5 


2.0 


GX 5-1 


9.0 


1.05 


166.0 


33.2 


110.0 


CygX-2 


8.0 


1.06 


11.5 


8.5 


7.3 


XB 1916-053 


10.0 


1.62 


0.66 


0.87 


0.08 


XB 1323-619 


10.0 


1.77 


0.26 


0.45 


0.02 


X 1624-490 


15.0 


1.31 


25.5 


8.3 


6.8 


XBT 0748-676 


10.0 


1.99 


0.94 


0.7 


0.05 



Thus, the large variation in emitting area is responsible 
for the large change in blackbody level. At low luminosity, 
the emission is from an equatorial strip on the star, which 
increases in height till at the higher luminosities, all of the 
star is emitting. 




1 10 
Rbb C^i") 



Fig. 3. Variation of blackbody temperature fcTBB and ra- 
dius i?BB for the survey sources plus additional sources 
discussed in the text. 90% confidence limits are shown 



2.6. Accretion Disk Properties 

We next express the blackbody emitting area in terms 
of the half-height h of the emitting strip assumed to be 
an equatorial strip of varying height. The area is 4 7r/ii?* 
(for a sphere intersected by two parallel planes) equal to 
47ri?|g. Thus, the half- height h is related to i?BB by 



^ = Rbb/R*- 



(1) 
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Fig. 4. Left: variation of the height h of the blackbody emission region on the surface of the neutron star with the 
equilibrium height i?eq of the radiatively supported inner accretion disk. The triangles show the blackbody radius 
i?BB instead of h for the two bright sources GX5-1 and GX13+1. Right: h of the blackbody emission region versus 
the more accurate accretion disk height H (see text) 



Values of h obtained in this way arc included in Table 5. 

Next, we compare h with the half-height of the in- 
ner radiatively-supported accretion disk. For a wide range 
of mass accretion rates, the inner disk is radiatively sup- 
ported with half-height H given by 



H 



3aT M 
87rmp c 



l_ 

r 



1/2. 



H, 



eq 



1- 

r 



1/2. 



(2) 



where ctt is the Thomson cross section, M is the mass 
accretion rate and nip is the proton mass (Frank, King & 
Rainc 1992). The disk height His zero at the surface of the 
star, but increases rapidly until it becomes independent of 
r having the equilibrium value H^q. Most of the increase 
takes place between r = 10 and r = 20 km, i.e. within 
10 km of the stellar surface. In weak sources, with L ~ 
5 X 10^^ erg s~^, radiative support is over a limited radial 
extent so that i?eq is never achieved. In bright sources, H^q 
is achieved, and radiative support continues typically to a 
radial distance of 400 km (for Lxot = 10'^* erg s~^). Values 
of Hgq are, of course, much larger than the height of the 
gas-dominated thin disk calculated using the solution of 
Shakura & Sunyacv (1976) which for luminosities between 
10^^ and 10^^ erg s~^, gives a half- height at a radius close 
to the neutron star (1.01 • i?,) varying between ~ 25-50 
m. 

As a first approximation, in Fig. 4 (left), we show the 
variation of h with H^q, which was calculated using the 
total luminosity of each source to give M via ixot = 
GM^,M / Rt:. It can be seen that there is agreement be- 
tween h and iJeq! however, the lower luminosity points 
fall below the line h = Heq. This would be expected, since 
for lower luminosities i?oq is never achieved since the disk 
is radiatively-supported over an insufficient radial extent. 



We calculate the actual maximum height of the disk in 
such cases as follows. The disk continues to be thick while 
the radiation pressure Pr » Pg, the gas pressure. Czerny 
& Elvis (1987) discuss the transition from a thin disk to a 
thick disk. From their Eqn. 19, it can be seen that the disk 
height becomes 90% of Heq for p^/ (pr +Pg) = 10. For a cal- 
culation of disk height accurate to ^^10%, we have derived 
the radius rio at which p^=10 pg, using the expressions for 
radiation pressure and gas pressure given by Shakura & 
Sunyaev (1976), to be no = 10.1 Mi"^/^^a2/2i ^16/21^^4/21 
km, where Mi is the mass of the neutron star in Solar 
masses, a is the viscosity parameter, Miq is in units of 
10^6 g s-i and / = [1 - (i?,/r)^/2]i/4 rj^j^jg equation can 
be solved numerically for rio for particular values of M. 
For example, using M values of 10^^ and 10^^ g equiv- 
alent to total luminosities of 2 x 10^'^ and 2 x 10"^® erg s~^, 
and assuming a = 0.1, rio is ~50 km and ~300 km respec- 
tively, but for 3 X 10^^ erg s~^, rio = 10.4 km only. The 
actual maximum height of the disk will be H{rio). In Fig 4 
(right) we show the variation of h with ff(rio), the more 
realistic calculation of disk height. It can be seen that 
there is remarkably good agreement between the height 
of the emission region and the height of the radiatively- 
supported disk, most points lying close to the line h = 
H. 

In the case of GX5-1, the agreement with h = H is 
poor, as h= 110 km and H = 43 km. However the black- 
body radius in this case is 33 km showing that the emis- 
sion is from a sphere three times larger than the star. 
There is a clear implication that the accretion flow is so 
strong that a spherical cloud of matter around the star has 
formed which is responsible for the blackbody emission. If 
we compare H with the height of this sphere, i.e. Rbb 
= 33 km, there is much better agreement between these 
quantities, again suggesting that the height of the disk de- 
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termines the height of the emission region. In Fig. 4 we 
plot triangles for GX5-1 and the other very bright source 
GX13+1 to show their expected positions a.t h = Rbb- 
Overall, there is good agreement between h and H in Fig. 
4 (right) over about 3 orders of magnitude on each axis 
which is strong evidence for the significance of the agree- 
ment. Appropriate errors are shown for all points derived 
from 90% confidence limits in spectra parameters. Error 
in source distance will affect both h and H; adopting 10 
km for the stellar radius will affect K a value of 12 km 
reducing heights by 20%. It is also clear that we may in- 
troduce error by using the luminosity in the band 1-30 
keV for calculating H, but this range was chosen to avoid 
extrapolating ASCA results too far above 10 keV. Using 
our work on XB 1916-053 from BeppoSAX (Church et al. 
1998a) we have estimated possible errors. For this source, 
with Eco = 80 keV, and the total luminosity would have 
been underestimated by 28% using the band 1-30 keV. 
For Eco = 20 keV, the error falls to 8% and for Eco 
= 10 keV, the error is 1%. Thus there can be a further 
source of error in H requiring points to be moved to the 
right by up to 30%. Clearly, broadband spectral analysis 
is desirable for all sources in the survey. The implications 
of Fig. 4 will be discussed in the next section. 

3. Discussion 

We have shown that the two-component blackbody plus 
Comptonization model gives a good description of the 
ASCA spectra of the Atoll and Z-track sources analysed 
in this survey. We thus now propose that this model sug- 
gested by us as a unifying model for the dipping sources 
(Church & Balucinska-Church 1995) gives a good descrip- 
tion of LMXB in general. Fitting results for the disk black- 
body plus Comptonization model were in 8 out of 13 cases 
unphysical having values of inner accretion disk radius 
substantially smaller than the neutron star radius making 
it unlikely that the blackbody emission originates in the 
disk. Our previous results for the dipping sources in par- 
ticular (see Sect. 1), prove that the Comptonizing region is 
extended, indicating a flat, vertically shallow ADC above 
the inner accretion disk. Thus we do not expect signifi- 
cant Comptonization of blackbody emission from the neu- 
tron star in such an ADC since it is geometrically unlikely 
that blackbody photons will pass through this region. It is 
similarly unlikely that substantial Comptonization of the 
neutron star emission takes place in the inner radiatively- 
supported disk as this would require long path lengths of 
photons through the disk. 

Next, we have shown that use of the two-component 
model reveals a broad dependence of the blackbody lumi- 
nosity on the total luminosity not previously known for 
LMXB in general, the blackbody luminosity falling more 
rapidly than the total luminosity as the mass accretion 
rate decreases. The major change in blackbody emission 
has been shown to result primarily from changes in the 



emitting area, not the temperature. Assuming that the 
emission originates in an equatorial belt on the neutron 
star having radius --^lO km, we have shown that its height 
agrees well with the height of the radiatively-supported 
inner accretion disk. This agreement extends to the less lu- 
minous sources where the radiatively-supported disk does 
not achieve its equilibrium height, and there is good agree- 
ment between h and H calculated at the position where 
Pt = 10 Pg. 

Possible explanations of the result that h = H are of 
two types; those in which the accretion disk height directly 
determines H, and those in which it does not. The direct 
mechanism requires radial (advcctive) flow of material be- 
tween the inner disk edge and the star. Although there has 
been a very large amount of theoretical work on advective 
flow in black hole binaries (e.g. Abramowicz ct al. 1996, 
and references therein), little work has been carried out 
for LMXB. It is not known whether the sonic point lies 
within the inner disk or within the star. However, Popham 
& Sunyaev (2000) show that the radial velocity increases 
by two orders of magnitude in the inner disk, but it is not 
known whether accretion flow can cross the gap between 
the inner disk edge and the star. Resolving this question 
will require detailed two- or three-dimensional hydrody- 
namic modelling. 

Inogamov & Sunyaev (1999) have recently shown that 
accretion flow meeting the surface of the neutron star in 
the equatorial plane will spread over the surface of the 
star producing regions of enhanced X-ray brightness of 
vertical extent increasing with luminosity, as is the case 
in our observational result. The mechanism proposed is 
independent of the accretion disk height and so the h = 
H agreement would result from H being a measure of L. 
Preliminary work shows encouraging agreement between 
the emitting height on the neutron star predicted by In- 
ogamov and Sunyaev. A detailed comparison between the 
present results and this theory will be made in a further 
paper. 

The above results have a bearing on alternative models 
for LMXB. Our results are very inconsistent with models 
in which the blackbody originates in the accretion disk, 
and the agreement of h = His positive evidence for neu- 
tron star blackbody emission. Additionally, we have used 
the timescales for dip ingress and egress in several obser- 
vations of dipping sources (Church 2000) to obtain the 
radius of the extended Comptonizing ADC region: tadC: 
typically 5 x 10^ cm. These measurements are supported 
by the fact that the Comptonized emission component is 
gradually covered during dipping. Thus the ADC is 5,000 
times larger than the neutron star, and this rules out mod- 
els in which Comptonization is supposed to take place in 
the neighbourhood of the neutron star. 

We clearly need to improve our knowledge of the emis- 
sion parameters of some of the sources studied by obtain- 
ing broadband spectra extending to ~100 kcV, allowing 
the cut-off energy to be measured, which has not in general 
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been possible in the present work. Further investigation 
of both very bright sources and very weak sources will 
also be of interest. A complete description of the black- 
body and Comptonizcd emission in LMXB must be more 
complicated than given here. For example, during flaring, 
major changes take place on the neutron star and inner 
disk and this is investigated in work currently taking place 
(Baluciiiska-Church et al. 2001). 

Acknowledgements. MJC and MBC thank the Royal Society 
and the British Council for financial support, and Andrew 
King, Wlodek Kluzniak, Nail Inogamov and Rashid Sunyaev 
for helpful discussions. 

References 

Abramowicz M.A., Chen X.-M., Granath M., Lasota J.-R, 

1996, ApJ 471, 762 
Asai K., Dotani T., Nagase F., Mitsuda K., 2000, ApJS 131, 

571 

Balucinska-Church M., Belloni T., Church M.J., Hasinger G., 

1995 A&A 302, L5 
Balucinska-Church M., Church M.J., Oosterbroek T. et al., 

1999 A&A 349, 495 
Balucinska-Church M., Humphrey P. J., Church M.J, Parmar 

A.N., 2000, A&A 360, 583 
Balucinska-Church M., Barnard R., Church M.J, Smale A. P., 

2001, A&A, submitted 
Barret D., Olive J.F., Boirin L. et al., 2000, ApJ 533, 329 
Christian D.J., Swank J. H., 1997, ApJ Suppl 109, 177 
Courvoisier T. L-J., Parmar A.N., Peacock A., Pakull M., 1988, 

ApJ 309, 265 

Church M.J., 2000, Proc. of 33rd Scientific Assembly of 
COSPAR, Warsaw, July 2000, Ed. Czerny B., Advances 
in Space Research, in press 
Church M.J., Baluciiiska-Church M., 1993, MNRAS 260, 59 
Church M.J., Balucihska-Church M., 1995, A&A, 300, 441 
Church M.J., Dotani T., Balucinska-Church M. et al., 1997, 
ApJ, 491, 388 

Church M.J., Balucinska-Church M., Dotani T., Asai K., 

1998a, ApJ, 504, 516 
Church M.J., Parmar A.N., BalucinskarChurch M. et al., 

1998b, A&A, 338, 556 
Czerny B., Elvis M.S., 1987, ApJ 321, 305 
Czerny B., Czerny M., Grindlay J.E., 1986, ApJ 311, 241 
Frank J., King A.R., Lasota J.-R, 1987, A&A 178, 137 
Prank J., King A.R., Raine D., 1992, "Accretion Power in As- 
trophysics" , Cambridge University Press 
Pukazawa Y. Ishida M., Ebisawa K., 1997, ASCA News, no. 4 
Guainazzi M., Parmar A.N., Segreto A., et al., 1998, A&A 339, 
802 

Hasinger G., van der Klis, 1989, A&A 225, 79 
Hasinger G., Priehorsky W.C., Middleditch J., 1989, ApJ 337, 
843 

Hasinger G., van der Klis M., Ebisawa K., Dotani T., Mitsuda 

K., 1990, A&A 235, 131 
Inogamov N.A., Sunyaev R.A, 1999, Astron Lett 25, 269 
in't Zand J.J.M., Heise J., Kuulkers E. et al., 1999, A&A 347, 

891 

Kubota A., Tanaka Y., Makishima K. et al., 2000, PASJ 50, 
667 



Merloni A., Fabian A.C., Ross R.R., 2000, MNRAS 313, 193 
Mitsuda K., Inoue H., Nakamura N., Tanaka Y., 1989, PASJ 
41, 97 

Orosz J. A., Kuulkers E., 1999, MNRAS 305, 132 

Parmar A.N., White N.E., Giommi P., Gottwald M., 1986, ApJ 

308, 199 

Parmar A.N., Gottwald M., van der Klis M., van Paradijs J., 

1989, ApJ 338, 1024 
Schulz N.S., 1999, ApJ 511, 304 

Schulz. N.S., Hasinger G., Triimpcr J., 1989, A&A 225, 48 

Shakura N.I., Sunyaev R.A., 1976, MNRAS 175, 613 

Smale A.P., Mukai K., Williams O.R., Jones M.H., Corbett 

R.H.D., 1992, ApJ 400, 330 
Smale A.P., Done C, Mushotszky R.F., et al, 1993, ApJ 410, 

796 

Sunyaev R.A., Shakura N.I., 1986, SvAL 12, 117 
Tanaka Y, Inoue H., Holt S.S., 1994, PASP 46, L37 
Titarchuk L., 1994, ApJ 434, 313 
White N.E., Swank J.H., 1982, ApJ 253, L61 
White N.E., Peacock A., Taylor B.C., 1985, ApJ 296, 475 
White N.E., Stella L., Parmar A.N., 1988, ApJ 324, 363 
van der Klis M., 1995, in Lewin W.H.G., van Paradijs J., van 
den Heuvel E. P. J., Eds., X-ray Binaries. Cambridge Uni- 
versity Press, Cambridge 



